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This work presents tables to assist in the evaluation of pinhole collection effects on spacecraft. These 
tables summarize results of a computer model which tracks particle trajectories through a simplified electric 
field in the plasma sheath. A technique Is proposed to account for plasma sheath efFects in the application 
of these results and scaling rules are proposed to apply the calculations to specific situations. This model is 
compared to ion current measurements obtained by another worker, and the agreement is very good. 


INTRODUCTION 

Ionospheric plasma interactions can result in 
spacecraft structure potentials being substantially 
more negative than the ambient environment 1 , 
especially on high power spacecraft supported by 
high voltage solar arrays. Among other things 
this can result in sputtering of spacecraft surfaces 
by the collected ions if the potential difference 
is negative enough. This sputtering may become 
a concern in three ways: it may be a source of 
contamination for surfaces closer to the plasma 
potential; it may cause failure of grounded metal 
foils; and as the effective ion collection area grows 
structure currents and power losses will grow. 

Sites of special concern are grounded metal 
surfaces covered with insulating dielectric films ex- 
posed to the ionospheric plasma. Even though in 
general these films will be protective, defects in 
these films caused by the manufacturing process, 
micrometeoroid and debris hits, or dielectric break- 
down, will result in the focusing of the collected 
ions onto local spots of the metal. The enhanced 
ion collection will be significantly more than that 
estimated from only the ion flux to the surface 
and the area of the defect, possibly resulting in 
enhanced sputter rates. 

In order to calculate the ion current to the 
surface a focus factor can be estimated. This 
focus factor is a characteristic of the geometry 
and conditions of and near the defect. Once it is 
evaluated, the current to the defect is given by the 
product of the ion flux and the defect area with 
the focus factor. 

This paper is declared a work of the U S. Government 
and is not subject to copyright protection in the United 
States. 
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Work on problems with similar electric field 
geometries has been performed by various workers. 
These problems include electron “snapover” near 
pinholes 2 ” 6 , and focusing onto dielectric surfaces 7 . 
These studies include physical processes that are 
not necessary for the ion collection analysis, though 
generally the final results are similar. 

The work presented here is an extension of an 
effort 8 to develop a simple quantitative model to es- 
timate the ion current collection to holes in space- 
craft insulating surfaces in ionospheric plasma. 
This work incorporates a technique to estimate the 
size of the plasma sheath, rather than treating the 
sheath as a free parameter. Also suggestions are 
made as to how the parameters can be scaled to 
permit application of a relatively small data set 
to a wide range of physical situations. Results 
of focus factor estimates are tabulated, permitting 
interpolation to estimate current collection. In ad- 
dition, calculated results are compared to some 
experimental observations. 

MODEL 

Original Model 

The approach taken in this model is to solve 
Laplace’s equation for a pinhole geometry in a 
region representative of a plasma sheath, then track 
ion trajectories through the resulting potential 
field. The fields and trajectories are determined in 
two dimensions with cylindrical symmetry assumed 
about the axis of the hole, as shown in Fig. 1. 
However velocities are tracked in three dimensions 
to permit accommodation of angular momentum 
considerations. The focus factor for a particular 
incident velocity vector is the ratio of the area of 
initial positions from which ions will hit in the hole 
to the area of the hole. Ion current to the hole can 
then be calculated from the product of the hole 
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area, the focus factor and the incident ion current 
density for the sheath edge. 

The calculation region used is representative 
of the plasma sheath. The electric fields near the 
hole are calculated by solving Laplace’s equation 
rather than Poisson’s equation. This is appropriate 
when the bias is large compared to the plasma 
temperatures, and the density of the collected ions 
does not get too high. This is usually true since 
the ions accelerate as they leave the plasma sheath 
edge. Also, while ion trajectories may be affected 
by increased ion density due to the focusing, in 
practice the ions tend to hit near the center of the 
conductor, and deflecting the trajectories outward 
slightly will not effect the observed ion current. 

The boundary conditions are simplified by as- 
suming the upper boundary of the rectangular 
calculation space is at zero potential. This rep- 
resents the upper boundary of the plasma sheath, 
and it is assumed to be parallel to the surface. 
The lower boundary is divided into two regions, a 
central hole region at the biased potential, and a 
dielectric region at zero potential. The bias po- 
tential represents the electric potential difference 
between the ambient plasma and the spacecraft 
ground. Since dielectrics generally charge to a 
small negative potential, comparable to the elec- 
tron temperature, and the hole region is metal, 
this assumption is reasonable as long as the bias is 
much larger than the potential difference attained 
between the dielectric and the plasma. The cen- 
tral and outer boundaries are set so that the radial 
electric field is zero. At the center this is true by 
symmetry. It is also true at the outer boundary if 
it is far enough away. 

This simplified field approach ignores the elec- 
tric field structure at the hole dielectric edge and 
the depth of the hole, which should have only short 
range effects. In this case short range means that if 
the particle gets close enough to experience them, 
it will probably hit within the hole. But the is- 
sue may be important for very small defects where 
the insulator conductivity may result in an effec- 
tive hole size substantially larger than the defect. 
It also ignores the electric field between the di- 
electric and the plasma due to the thermal plasma 
sheath. However, this should be negligible, and the 
assumption reduces the dimension of the model s 
parameter space. 

Secondary electron emission effects are ne- 
glected. ’Snapover’ on the insulation is not im- 
portant in this case since the electric field on the 
insulators prevents the secondaries from escaping. 



Figure 1: Calculation Region 


On the conductor, secondary yields may be on 
the order of 10%, but this is small compared to 
uncertainties in the plasma parameters and may 
be neglected for engineering level calculations, or 
added explicitly on a case by case basis. 

The particle velocity is tracked in three di- 
mensions. The potential field around the hole has 
cylindrical symmetry so the position angular co- 
ordinate of the ion trajectory is not important. 
However the angular velocity component plays a 
significant role since this component is needed to 
conserve angular momentum about the axis of 
symmetry. Angular effects are accounted for by 
calculating the displacement over the time inter- 
val, then rotating the coordinate system to find 
the new values for r, z y tv, v z and v$. 

The focus factor is the ratio of the area where 
incident particles will hit inside the hole, to the 
area of the hole. For cylindrical symmetry this is 
the ratio, r}/r\ y where is the maximum incident 
radius that ions are collected from, and r* is the 
radius of the hole. At large angles of incidence, this 
expression may have to be corrected to account for 
the fact that there can be a minimum radius, in- 
side which incident ions are not collected. This 
definition is convenient for determining ion cur- 
rents, and sputter rates as contamination sources. 
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Figure 2: Illustration of incident velocity varia- 
tion for the different sectors of the calculation. 

However, this is an underestimate for local dam- 
age rates since almost all the collected particles hit 
inside the radius, r*/ 2. 

Even though the electric field in this prob- 
lem has cylindrical symmetry, this symmetry is 
destroyed by an incoming flux of particles with a 
non-zero angle of incidence, In this case the fo- 
cus factor depends on the angle coordinate, 9. In 
order to take into account the 9 dependence, the 
360 degree range for values of 9 are divided into 
a number of sectors as illustrated in Fig. 2, each 
with a representative value for 9 , and a represen- 
tative value of the focus factor is found for that 
sector. An even number of sectors is used, includ- 
ing ones centered on 0 and 180 degrees. In practice 
45 degree sectors are used. However, because of 
the symmetry of the problem only half the other 
sectors need to be calculated, either 0 < 9 < 180 
degrees or 180 <9 < 360 degrees. Those sectors at 
zero and 180 degrees are given a weighting of one 
while the others are given a weighting of two to 
account for the missing but opposite sectors. The 
net focus factor is the average of the sector focus 
factors. 

The effect of an ion temperature can be ad- 
dressed by adding incremental velocities repre- 
sentative of a three-dimensional Maxwell velocity 
distribution to an incident velocity. To calcu- 
late a focus factor for a given incident angle and 
temperature, the three dimensional initial velocity 


distribution is broken into a cartesian grid. A rep- 
resentative velocity vector is chosen at the center 
of each grid element and a focus factor is calcu- 
lated for each grid element. Since only half the 
elements need be calculated because those with a 
velocity element in the -i-0 have the same focus 
factor as the corresponding element in the —9 di- 
rection, division of the three dimensional velocity 
space into approximately 3000 elements (18x18x9) 
yields a reasonable simulation. Each grid element 
is assigned a composite velocity vector represen- 
tative of its position and weighted first by the 
probability of the particle having a velocity near 
the tested value, and secondly by which 9 sector 
the computation is being performed in. 


Plasma Sheath 

Since in the present model the plasma sheath 
region is equated with the calculation space, it 
is critical to be able to estimate the thickness of 
the sheath. This is accomplished by comparing 
the Laplacian solution with no upper boundary 
to the charge density of the influx of ions from 
the plasma. The location of the sheath edge, h s 
can be estimated by determining where the second 
partial derivative of the potential, <£, with respect 
to height above the hole is equal to the ion charge 
density flowing into the sheath. This condition is 
expressed by, 


d 2 <j> _ p(h,) 

dh* h , co • 


This is not a rigorous technique but finds the place 
where the geometric distortion of the field is com- 
parable to that due to the space charge. If the 
magnitude of each of the two terms of the Lapla- 
cian operator for cylindrical symmetry is much 
larger than the charge density term, the poten- 
tial field is primarily influenced by the boundary 
conditions, and Laplace’s equation is appropriate. 
But when the terms get as small as the charge 
density term, Poisson’s equation clearly should be 
used. We avoid the problem by assuming this 
condition indicates the edge of the plasma sheath. 
This condition is used to estimate the sheath edge 
location. 
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Figure 3: Comparison between an analytical ex- 
pression, V a = (1 — h/yffTT?), for the electric 
potential above a pinhole and a calculated finite 
difference solution, V C) to Laplace's equation. 


In order to do this comparison, an analytical 
expression was found for the potential along the 
axis above the hole. A technique was used which 
has occasionally, but not reliably, found analyti- 
cal solutions to Laplace’s equation. The potential 
at a point is assumed to be the average of defined 
surface potentials weighted by the solid angle cov- 
ered. The result is then checked against Laplace’s 
equation to see if it is correct. The average an- 
gular potential, Va, seen at a height, h , above the 
center (r=0) of the biased hole is obtained from, 

^arctanir^jh) 

V a (h) = / <t> (ft=0, r) sin (0) dd 

where <£(/i=0, r) = V 0 for r < r h and <£(/i=0,r) = 0 
for r > tv Since this expression was not found 
rigorously, it is compared to a solution found 
by finite differencing in Fig. 3. The maximum 
difference between the two is about 0.001 times the 
bias voltage within 10 hole radii on the conductor. 
It can be seen that the agreement is excellent until 
the finite difference solution is forced to zero at the 
upper boundary of the calculation space. 

The value obtained above should be compared 
with the Debye length for the plasma. The Debye 
length gives an idea of the plasma sheath associated 


with the surface, and if the surface’s thermal 
sheath is thicker than the sheath associated with 
the hole, the assumptions made in this model are 
not appropriate. 


Scaling Rules 

The pinhole ion collection problem has six pa- 
rameters; hole radius and bias, plasma speed and 
temperature, the angle that the ions’ velocity vec- 
tor makes with the axis of symmetry, and the 
plasma sheath thickness. However, it is apparent 
that a number of interesting parameters are inter- 
related so that results are reported in terms of only 
four parameters: normalized plasma ram energy 
and incident angle, normalized ion temperature, 
and normalized sheath thickness. These simplifica- 
tions are reasonable if two issues are demonstrated, 
(1) the electric field depends only on the grid cho- 
sen and can be scaled with the magnitude of the 
applied bias and the hole radius, and (2) the tra- 
jectory of the particle within the grid depends only 
on the scaled electric field, the incident kinetic en- 
ergy to bias ratio, and the velocity direction. The 
following discussion, while not rigorous, suggests 
that this is indeed the case. 

The electric field in the calculation space is 
obtained by setting certain boundary conditions, 
and satisfying Laplace’s equation. First consider 
only the applied bias. If V 2 <ji(r, z) = 0 then 
V 2 Vo <j>{r>z) = 0 is also satisfied, where Vq is a 
constant adjusted to match the voltage difference. 
The electric field is given by V 0 V^*(r, z) where 

z) is the potential field that matches the 
boundary conditions when the bias is one. 

There are three distance scales in the defini- 
tion of the problem, the radius of the hole, the 
height of the calculation space, i.e. the plasma 
sheath thickness, and the width of the calculation 
space. Ideally the width of the calculation space 
would be infinite, and in practice it should be large 
enough that it does not effect the calculation sig- 
nificantly. This leaves only two distance scales, the 
sheath thickness in the z direction and the hole 
radius in the r direction. 

In practice the potential field is found for po- 
sitions on a grid using a finite difference equation. 
For cylindrical symmetry, this comes from, 
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At the grid position (t'Ar, jAz) the finite difference 
equation becomes, 




(t + <fc+l,j 4- («' — I) 


2 (a 2 + 1) 


I ' O 

a 1 


where a = Az/Ar . Note that the grid separation 
drops out. As long as the grid is cartesian the 
potential field depends only on the boundary con- 
ditions and the relative grid spacing. So as long as 
the width of the grid range is enough to be effec- 
tively infinite, the potential field, which depends 
only on the boundary conditions, depends only on 
the plasma sheath thickness height and hole ra- 
dius ratio, for the relative grid size, and the bias 
potential. 

It is the intention here to show that the 
electric field can be scaled to the grid size or the 
hole radius. An electric field is obtained from the 
potential grid using equations like, 


E r ( iArJAz ) = - Vo 


2Ar 


Note that for a hole radius of r* = (/ 4- |)Ar, a 
’grid’ electric field, E*(i, j), relative to the hole 
size can be defined as, 


™ .x _ r h E r (iArJAz) 
E - = Vo 



Similarly E*(iJ) can be obtained. The scaled 
electric field, (V 0 /r h )E*(iJ)> is then used in the 
equations of motion for the electric field. 

Now the motion within the grid is examined 
assuming that the electric field can be scaled by 
E(r,h) = V 0 E*(r/r hi h/r h )/r h . The model uses 
the two equations of motion: 


small compared to E. But this time interval can 
be chosen so that 6t = 6x/|vb|, where 6x is a 
distance interval along the particle path. Substi- 
tuting the scaled electric field above, the kinetic 
energy, K = ^m|tT| 2 , and scaling to the hole ra- 
dius distance scale, n, produces the following two 
equations of motion: 

v(6t) _ |iTo| ( vo , qV66x ( r_ h\\ 

r h ~ r* \\v 0 \ 2 Kr h \r h ' r h J ) ’ 

r(6t) _ r 0 vpSx gVpSx 2 / _r_ h\ 

r h |t? 0 |r* AKr\ \r h ' r h ) ' 

The scaled trajectory, obtained from the position 
equation depends only on the direction of the 
velocity vector, vo/|t7o|> the energy ratio, qVp/K , 
and the scaled electric field E*(r/rh , h/ru). The 
accuracy of the equation improves as the scaled 
interval Sx/r^ gets smaller. 



Figure 4: Focus Factor for a 3 mm radius hole as 
a function of angle under Low Earth Orbit (LEO) 
conditions as described in the text. 


= v 0 + a6t , 
r(6t) = ro + vodt "V 

where a is the acceleration is given by qE(r,h)/m 
and 6t is a time interval to increment the cal- 
culation. In the limit that 6t goes to zero, this 
will converge on the correct trajectory, jand will 
give good results when (f(6t) — f 0 ) ■ VE is very 


RESULTS 

Normalized results of the computer model are 
tabulated in the Appendix. The tables show the 
focus factor as a function of the four variables; 
S', ©', E\ and S' is the plasma sheath 
thickness divided by the hole radius. 0', is the 
plasma ion temperature in eV divided by the 
absolute value of the hole bias in Volts. E\ is the 
directed incident kinetic energy of the ions in eV 
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divided by the absolute value of the hole bias in 
volts. Finally, is the angle between the hole’s 
axis of symmetry and the incident velocity vector, 
ignoring temperature. ($ is a characteristic of the 
population of ions, whereas <p is a characteristic of 
individual ions.) 

Focus factors for low energies and high angles 
of incidence should be treated with caution. For 
low energies focus factors tend to be high. It 
appears that all particles entering the grid with 
zero velocity are collected by the hole, giving an 
infinite focus factor contribution. In this case 
the electric field due to the thermal sheath must 
be included to bound the focus factor. In the 
tables presented here focus factors in this case 
are bounded by the size of the calculation space 
instead. 



square root of K/Ti 


Figure 5; Effect of a thermal distribution on aver- 
age velocity of incident ions as the ratio, Kzo/qTi 
varies. The velocity is normalized by the factor 
y/m/2qTi. 


At large angles of incidence there is also a 
possibility that the focus factors are artificially 
high. In this case particles which hit the hole 
come from far away and receive an area weighting 
proportional to 1/r. However, nearby particles 
with a slightly different vg component may not be 
collected due to higher angular momentum. In this 
case the coarse sector divisions used may result in 
undue weight given to the large focus factor. 

These tables can be interpolated for specific 
situations. As a example, Fig. 4 shows the focus 


factor as a function of incident angle, for a 
3 mm radius hole biased to -80 V in a typical 
low earth orbit plasma. This oxygen ion plasma 
has a temperature of 0.1 eV, a density of 1 x 
10 6 cm~ 3 , and an incident kinetic energy of 4.5 eV. 
In this case the plasma sheath is found to be 
19 mm thick so S' = 6.4, ©' = 0.0013, and E f — 
0.056. A polynomial interpolation scheme 9 was 
used to interpolate the ram energy and angle 
of incidence, and a cubic spline 9 was used to 
interpolate the temperature and sheath thickness 
parameters. The ability of the of the interpolation 
scheme to reproduce calculated results not included 
in the tables is within 2%. 

The focus factors from Fig. 4 show a slight 
decrease as the incident angle increases from zero, 
then a pronounced increase at higher angles of 
incidence. The initial decrease with angle may be 
due to the added angular momentum preventing 
particles from hitting the hole. At high angles 
of incidence the particles are coming from farther 
away. The increased area weight results in an 
enhanced focus factor. 



Figure 6: Ion currents to a 3 mm radius hole as a 
function of angle under LEO conditions, illustrat- 
ing effects of the ion focusing shown in Fig. 4. 

The focus factor, /, is converted into a current 
per hole, t, using f = n t q{v z )}A , where A is the 
area of the conductive hole, n e is the plasma 
density, q is the charge and (v z ) is the average ion 
velocity normal to the surface obtained from: 
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Figure 7: Comparison between Estimated Cur- 
rents and Measured pinhole ion currents. The 
upper set of curves are obtained for a radius of 
2.4 mm and the lower set is for a radius of 1.4 mm. 
The plasma conditions are described in the text. 


or, 



were v z o is the initial non-thermal velocity normal 
to the surface, K z p is the kinetic energy due to that 
motion and qTi is the ion thermal energy when Ti is 
the ion temper ature in e V. Fig. 5 shows that (v z ) 
converges to \J2 K z o/m very quickly. Assuming 
{ v z ) = y/2K z0 /m at y/K z0 /qTi = 1 is low by only 
2.5%. Figure 6 includes the current densities to 
give currents to the 3 mm radius hole. In this 
case the current stays near 0.3 /iA until very high 
angles of incidence where it drops off. 

COMPARISON TO EXPERIMENT 

The model has been compared with recent 
experiment measurements and is reasonably con- 
sistent with those results. However, the experiment 
tests only normal incidence to the hole, and does 
not verify the angular dependence. 

Vaughn 10 has measured ion currents as a func- 
tion of bias voltage to two differently sized holes. 
He measured the current to a 2.7 mm radius hole 


and a 1.4 mm radius hole biased at from 0 V to 
-230 V. The Argon plasma he used had an elec- 
tron temperature of 1.2 eV and an ion energy of 
2.0 eV. However, the ions were thought to be cold. 
The number density differed between the two sets 
of data. For the 2.7 mm hole the density was 
2 x 10 12 m“ 3 giving an estimated current density 
of 1 mA/m 2 , and for the 1.4 mm hole the den- 
sity was 3 x 10 12 m~ 3 giving a current density 
of 1.5 mA/m 2 . Figure 7 compares the measured 
currents to those obtained from the calculated fo- 
cus factor and current density. The estimated 
currents are shown for two cases, one with cold 
ions, Ti = 0 eV, and the other for warmer ions, 
Ti = 0.2 eV. Even for the cold case agreement is 
within 30% of the measured values over the -50 V 
to -200 V range, within the uncertainty of the 
plasma parameters. However, assuming warmer 
ions improves the agreement slightly. 

CONCLUSION 

The objective of this work is to produce a 
model which can be used to predict ion currents 
to insulator defects under ionospheric conditions. 
These currents are needed to estimate either space- 
craft floating potentials or material sputter rates. 

The simple approximations used in the devel- 
opment of the computer model described here are 
adequate to estimate ion currents to conductive 
defects on insulation. While the original intention 
for model development was to produce a model 
easily run on small computers, the addition of 
thermal distributions requires long computation 
times. However, the model has been reduced to 
depend on four parameters which can be tabu- 
lated. The Appendix contains tables which span a 
reasonable space for small holes in ionospheric con- 
ditions. These tables predict ion currents that are 
consistent with those measured during an initial 
test by other workers. 

This paper only addresses part of the issue, 
the part dealing with individual defects geometries. 
In order to estimate currents and sputter rates 
of spacecraft surfaces, estimates of the material 
specific defect populations as well as estimates of 
the sputter yield are also needed. 
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APPENDIX 


The 28 tables show the focus factor as a func- 
tion of the four variables; S', ©J, E-, and $. S' 
is the plasma sheath thickness divided by the hole 
radius and ranges from 3.6 to 11.2. ©J is the plas- 
ma ion temperature in eV divided by the absolute 
value of the hole bias in volts and ranges from 0 


to 0.2. E'f is the directed incident kinetic energy of 
the ions in eV divided by the absolute value of the 
hole bias and ranges from 0.01 to 0.10. Finally, $ is 
the angle between the hole’s axis of symmetry and 
the incident velocity vector, ignoring temperature, 
and ranges from 0 to 75 degrees. 


Table 1. $ = 0.0 E,' = 0 01 


S' ©' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

13.3 

13.9 

14.8 

15.8 

15.3 

12.7 

12.4 

4.1 

15.1 

15.8 

16.9 

17.9 

17.2 

14.1 

13.8 

4.9 

17.9 

18.8 

20.2 

21.4 

20.3 

16.4 

16.0 

6.4 

22.1 

23.1 

24.8 

26.5 

24.6 

19.5 

19.0 

11.2 

31.9 

34.2 

36.8 

37.2 

34.8 

25.2 

23.8 


Table 3. 

$ = 

0.0 e; = o.o3 


s' e; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

8.0 

7.8 

8.0 

8.6 

9.3 

9.4 

9.8 

4.1 

8.7 

8.6 

8.8 

9.5 

10.3 

10.3 

10.8 

4.9 

9.4 

9.9 

10.1 

10.9 

11.7 

11.8 

12.3 

6.4 

10.9 

11.4 

11.7 

12.7 

13.6 

13.5 

14.3 

11.2 

14.2 

14.7 

15.0 

16.2 

16.9 

17.3 

17.9 



Table 2. 

$ = 

0.0 E 

= o.o2 


S' ©' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

9.4 

9.7 

10.0 

11.0 

11.5 

10.7 

11.0 

4.1 

10.9 

10.8 

11.2 

12.3 

12.8 

11.9 

12.1 

4.9 

12.5 

12.6 

13.0 

14.4 

14.8 

13.6 

13.9 

6.4 

14.2 

14.9 

15.5 

17.1 

17.5 

15.8 

16.3 

11.2 

20.0 

20.2 

20.9 

22.9 

21.8 

19.0 

19.7 


Table 4. 

4* = 

0.0 L 

;; = o.io 


S' ©' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

4.0 

4.0 

4.0 

4.1 

4.2 

4.4 

4.6 

4.1 

4.0 

4.2 

4.2 

4.3 

4.5 

4.7 

4.9 

4.9 

4.5 

4.6 

4.6 

4.7 

4.9 

5.1 

5.4 

6.4 

5.0 

5.0 

5.0 

5.1 

5.3 

5.6 

5.8 

11.2 

5.5 

5.6 

5.6 

5.7 

6.0 

6.1 

6.5 
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Table 5. 

$ = 

5.0 E 

= 0.01 



Table 6. 

$ = 

5.0 E 

= 0.02 


5' e'i 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

5' e ; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

13.1 

13.8 

14.6 

15.2 

14.6 

12.1 

11.8 

3.6 

9.6 

9.7 

9.9 

10.7 

10.9 

10.1 

10.3 

4.1 

14.9 

15.6 

16.5 

17.2 

16.4 

13.5 

13.2 

4.1 

10.4 

10.7 

11.0 

11.9 

12.1 

11.1 

11.4 

4.9 

18.0 

18.6 

19.6 

20.3 

19.3 

15.7 

15.5 

4.9 

12.5 

12.4 

12.7 

13.7 

13.8 

12.7 

13.1 

6.4 

21.8 

22.6 

23.8 

24.5 

23.3 

18.4 

18.5 

6.4 

14.6 

14.6 

14.9 

15.9 

15.9 

14.3 

14.8 

11.2 

32.4 

32.2 

33.4 

32.0 

27.0 

20.0 

17.9 

11.2 

19.6 

18.9 

18.9 

19.4 

17.6 

14.8 

14.2 


Table 7. 

$ = 

5.0 E 

' { = 0.03 



Table 8. 

$ = 

5.0 E 

' = 0.10 


s' e ; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

s' ©;• 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

7.7 

7.8 

7.9 

8.4 

8.8 

8.8 

9.1 

3.6 

3.9 

3.9 

3.9 

4.0 

4.1 

4.2 

4.3 

4.1 

8.4 

8.5 

8.7 

9.2 

9.7 

9.6 

9.9 

4.1 

4.2 

4.2 

4.2 

4.2 

4.3 

4.4 

4.5 

4.9 

9.5 

9.7 

9.8 

10.4 

10.9 

10.8 

11.1 

4.9 

4.5 

4.5 

4.5 

4.5 

4.5 

4.6 

4.7 

6.4 

11.2 

11.2 

11.3 

11.8 

12.2 

11.9 

12.1 

6.4 

5.0 

4.8 

4.7 

4.7 

4.7 

4.8 

4.8 

11.2 

14.3 

13.6 

13.5 

13.7 

13.1 

12.3 

11.5 

11.2 

5.7 

5.0 

4.7 

4.5 

4.5 

4.3 

4.2 


Table 0. 

$ = 

15.0 J 

B9i 



Table 10. $ = 




s' e; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

s' e[ 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

13.6 

14.0 

14.5 

14.5 

13.4 

11.1 

10.8 

3.6 

10.0 

9.8 

9.9 

10.3 

10.0 

9.1 

9.0 

4.1 

15.5 

15.8 

16.3 

16.2 

14.8 

12.2 

11.7 

4.1 

11.0 

10.8 

10.9 

11.2 

10.8 

9.8 

9.5 

4.9 

18.3 

18.5 

19.0 

18.6 

16.8 

13.6 

12.8 

4.9 

12.8 

12.4 

12.4 

12.6 

11.9 

10.5 

10.3 

6.4 

22.7 

22.1 

22.4 

22.2 

18.8 

14.6 

13.7 

6.4 

15.0 

14.3 

14.1 

13.9 

12.7 

11.0 

10.8 

11.2 

33.3 

29.0 

27.7 

22.8 

17.9 

13.2 

11.9 

11.2 

20.5 

16.6 

15.4 

13.5 

11.1 

9.1 

8.9 


Table 11. $ = 

15.0 

E '. = 0.03 



Table 12. $ = 

15.0 

El = 0.10 


s' e'i 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 _ 

s' e'i 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

7.9 

7.9 

7.9 

8.1 

8.1 

7.8 

7.1 

3.6 

4.1 

4.0 

3.9 

3.9 

3.8 

3.8 

3.7 

4.1 

8.8 

8.6 

8.6 

8.8 

8.7 

8.3 

7.4 

4.1 

4.3 

4.2 

4.1 

4.0 

3.9 

3.8 

3.8 

4.9 

9.9 

9.7 

9.6 

9.6 

9.3 

8.8 

7.7 

4.9 

4.8 

4.4 

4.3 

4.1 

3.9 

3.7 

3.7 

6.4 

11.9 

10.9 

10.7 

10.3 

9.7 

9.0 

7.7 

6.4 

5.2 

4.5 

4.3 

3.9 

3.6 

3.4 

3.3 

11.2 

15.0 

11.8 

10.8 

9.3 

7.9 

7.4 

6.2 

11.2 

5.8 

3.9 

3.5 

3.0 

2.6 

2.3 

2.3 


Table 13. $ = 


EBEH 



Table 14. 4> = 

wmm 

Ei=i 



s' ©; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

S ' ©' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

15.5 

15.2 

15.2 

14.1 

10.6 

9.8 

7.3 

3.6 

11.2 

10.6 

10.5 

10.1 

9.2 

8.0 

6.6 

4.1 

17.6 

16.9 

16.8 

15.4 

11.3 

10.4 

7.7 

4.1 

12.4 

11.6 

11.4 

10.9 

9.7 

8.4 

6.8 

4.9 

20.8 

19.4 

19.2 

17.2 

12.2 

11.2 

7.9 

4.9 

14.4 

13.1 

12.7 

11.7 

10.3 

8.8 

6.9 

6.4 

25.3 

22.7 

21.9 

19.2 

12.8 

12.0 

8.0 

6.4 

17.0 

14.6 

13.8 

12.2 

10.7 

9.1 

6.8 

11.2 

37.2 

26.4 

22.8 

17.0 

10.5 

9.7 

6.2 

11.2 

23.0 

15.1 

13.4 

10.5 

8.2 

7.6 

5.0 


Table 15. $ = 

30.0 

E\ = 

0.03 



Table 16. <t> = 

30.0 

% = 

0.10 


S' e'i 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

5' ©' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

9.1 

8.5 

8.4 

8.2 

7.7 

7.0 

6.1 

3.6 

4.7 

4.3 

4.1 

3.9 

3.6 

3.4 

3.3 

4.1 

10.2 

9.3 

9.0 

8.7 

8.0 

7.2 

6.2 

4.1 

4.9 

4.4 

4.2 

3.8 

3.5 

3.3 

3.2 

4.9 

11.5 

10.2 

9.8 

9.2 

8.3 

7.3 

6.3 

4.9 

5.3 

4.5 

4.2 

3.7 

3.4 

3.2 

3.0 

6.4 

12.9 

11.1 

10.3 

9.2 

8.1 

7.3 

6.1 

6.4 

5.4 

4.1 

4.0 

3.7 

3.3 

3.1 

2.9 

11.2 

15.5 

11.2 

10.2 

8.2 

6.5 

5.9 

4.9 

11.2 

7.2 

6.9 

6.3 

5.0 

3.9 

3.4 

3.1 
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Ta ble 17. $ = 45.0 = 0.01 Table 18. $ = 45.0 El = 0.02 


s' e; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

S' ©' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

19.3 

17.8 

16.8 

13.7 

10.4 

7.5 

6.9 

3.6 

13.5 

12.4 

12.2 

10.8 

8.8 

6.9 

6.2 

4.1 

21.4 

19.7 

18.8 

15.0 

11.1 

7.9 

7.1 

4.1 

14.9 

13.6 

13.3 

11.6 

9.2 

7.1 

6.5 

4.9 

25.1 

22.7 

21.7 

16.7 

12.0 

8.2 

7.5 

4.9 

17.6 

15.3 

14.8 

12.5 

9.8 

7.4 

6.9 

6.4 

30.6 

26.1 

24.6 

18.4 

13.0 

8.3 

7.8 

6.4 

20.4 

17.0 

15.9 

13.2 

8.7 

7.6 

7.4 

11.2 

45.6 

31.9 

25.4 

15.7 

11.2 

7.2 

6.0 

11.2 

26.9 

23.6 

19.5 

12.7 

9.2 

6.7 

6.1 


Table 19 

3> = 

45.0 

E\ = 0.03 



Table 20. 

$ = 

45.0 

E'; = 0.10 


5' e; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

s' e; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

10.9 

10.0 

9.8 

9.1 

7.8 

6.5 

5.2 

3.6 

5.6 

5.0 

4.7 

4.4 

4.2 

4.0 

3.6 

4.1 

12.3 

10.9 

10.6 

9.6 

8.1 

6.6 

5.3 

4.1 

5.7 

4.7 

4.7 

4.6 

4.4 

4.1 

3.7 

4.9 

13.9 

11.9 

11.4 

10.2 

8.5 

6.9 

5.4 

4.9 

5.6 

5.5 

5.6 

5.4 

4.9 

4.4 

4.0 

6.4 

15.9 

12.8 

12.3 

10.8 

8.9 

7.3 

5.5 

6.4 

7.6 

7.4 

7.3 

6.5 

5.5 

4.9 

4.4 

11.2 

23.4 

20.2 

16.9 

11.4 

8.3 

6.9 

5.1 

11.2 

13.8 

12.6 

10.8 

7.9 

5.9 

5.2 

4.4 


Table 21. $ = 

60.0 

E , < = 0.01 



Table 22. $ = 

60.0 

E '< = 0.02 


5' e; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

S' 0' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

25.8 

22.4 

18.7 

13.5 

9.0 

7.6 

6.2 

3.6 

18.1 

17.4 

16.1 

12.2 

8.6 

6.0 

5.3 

4.1 

28.5 

26.0 

21.4 

15.0 

9.5 

8.0 

6.3 

4.1 

20.6 

19.8 

18.3 

13.6 

9.1 

6.2 

5.4 

4.9 

34.4 

31.8 

25.8 

18.5 

10.4 

8.7 

7.1 

4.9 

25.0 

23.5 

21.8 

15.3 

10.0 

6.5 

5.9 

6.4 

43.4 

39.5 

31.0 

19.4 

11.3 

9.7 

7.8 

6.4 

30.6 

28.9 

25.8 

17.2 

11.0 

7.3 

6.6 

11.2 

75.1 

46.3 

29.6 

17.0 

10.4 

8.3 

5.8 

11.2 

60.5 

36.2 

24.9 

15.5 

8.9 

6.3 

5.5 


Table 23. $ = 


n 



Table 24. = 

60.0 

E . f = 0.10 


s' e; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

S' ©' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

15.4 

14.5 

13.9 

11.3 

7.6 

6.3 

4.5 

3.6 

7.8 

8.1 

8.3 

8.0 

6.6 

5.4 

4.5 

4.1 

17.1 

16.3 

15.6 

12.6 

8.0 

6.6 

4.7 

4.1 

9.5 

9.7 

9.8 

9.1 

7.2 

5.8 

4.8 

4.9 

20.2 

19.1 

18.5 

14.1 

8.7 

7.4 

5.0 

4.9 

12.2 

12.3 

12.0 

10.2 

7.7 

6.3 

5.3 

6.4 

24.4 

24.8 

22.5 

15.9 

9.8 

8.2 

5.6 

6.4 

16.2 

15.6 

14.1 

10.8 

8.0 

6.5 

5.3 

11.2 

51.0 

31.2 

21.7 

13.7 

8.9 

7.2 

5.1 

11.2 

29.7 

20.4 

15.5 

10.1 

7.6 

6.3 

5.0 


Table 25. 4> = 

75.0 

imm 



Table 26 

$ = 


e: = o.o 2 


s' e; 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

S' ©' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

40.2 

24.6 

17.8 

12.9 

7.5 

6.1 

5.3 

3.6 

34.4 

24.1 

18.7 

11.6 

8.2 

5.1 

4.4 

4.1 

55.5 

30.5 

21.1 

14.6 

8.0 

6.5 

5.6 

4.1 

48.1 

30.4 

22.7 

13.4 

8.9 

5.4 

4.6 

4.9 

76.6 

40.9 

26.9 

17.0 

8.8 

7.1 

6.1 

4.9 

64.1 

41.6 

29.7 

15.3 

10.1 

6.0 

5.1 

6.4 

108. 

55.8 

32.9 

20.2 

9.9 

7.7 

6.8 

6.4 

92.4 

51.6 

34.9 

17.1 

11.2 

6.9 

5.8 

11.2 

211. 

44.2 

26.4 

15.8 

8.1 

6.6 

5.2 

11.2 

160. 

28.6 

20.1 

12.8 

9.7 

5.8 

4.9 


Table 27. $ = 

75.0 

E \ = 0.03 



Table 28. $ = 

75.0 

Ei - = 0.10 


S' 0' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

S' ©' 

0.000 

0.001 

0.002 

0.005 

0.010 

0.015 

0.020 

3.6 

31.4 

23.8 

19.3 

10.9 

7.8 

6.0 

4.0 

3.6 

27.3 

21.7 

17.7 

11.3 

7.6 

5.7 

4.4 

4.1 

43.6 

30.4 

23.8 

12.7 

8.3 

6.5 

4.3 

4.1 

32.7 

26.0 

20.0 

11.9 

7.8 

5.9 

4.6 

4.9 

57.6 

42.1 

30.7 

14.6 

9.3 

7.5 

4.8 

4.9 

40.4 

29.4 

21.2 

12.1 

7.8 

6.1 

5.1 

6.4 

92.4 

49.8 

32.1 

16.3 

10.4 

8.6 

5.4 

6.4 

51.6 

29.9 

20.3 

10.9 

7.5 

6.7 

5.5 

11.2 

132. 

22.7 

16.6 

11.3 

7.6 

6.6 

4.5 

11.2 

84.4 

16.7 

10.8 

7.8 

6.1 

5.2 

4.3 
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